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Abstract. Neural field models with transmission delay may be cast as abstract delay 
differential equations (DDE). The theory of dual semigroups (also called sun-star calcu- 
lus) provides a natural framework for the analysis of a broad class of delay equations, 
among which DDE. In particular, it may be used advantageously for the investigation 
of stability and bifurcation of steady states. After introducing the neural field model in 
its basic functional analytic setting and discussing its spectral properties, we elaborate 
extensively an example and derive a characteristic equation. Under certain conditions 
the associated equilibrium may destabilise in a Hopf bifurcation. Furthermore, two Hopf 
curves may intersect in a double Hopf point in a two-dimensional parameter space. We 
provide general formulas for the corresponding critical normal form coefficients, evaluate 
these numerically and interpret the results. 



1. Introduction 

Spatial coarse graining of neural networks leads to so-called neural field models in which 
the average firing rates of underlying populations of neurons, as opposed to individual 
neuronal spikes, are considered. Such models have not changed substantially since the 
seminal work of Wilson and Cowan [53, 54], Amari [2] and Nunez [43]. Due to intrinsic 
delays of axons, synapses, and dendrites in the natural system, the role of delays in 
spatiotemporal dynamics of neural activity has received considerable attention [41, 34, 
32, 44, 52, 33, 31, 15, 14, 13]. Faugeras and coworkers investigated stability properties 
of stationary solutions using methods from functional analysis [49, 27, 50]. A first step 
towards Hopf bifurcation is made in [48], where Hopf bifurcation curves are computed. In 
[25] the principle of linearised stability and the Hopf bifurcation to periodic orbits were 
studied in the absence of delays. 

To set the stage, we have in mind p > 1 populations consisting of neurons that occupy 
fixed positions in a non-empty, bounded, connected, open region £1 C M n . For each 
i = 1, ... ,p let Vi(t, r) be the membrane potential at time i, averaged over those neurons 
in the ith population positioned at r £ ft. These potentials are assumed to evolve in the 
absence of time dependent external stimuli according to the system of integro-differential 
equations 

-^i(t,r) = -a i V i (t,r) + J2 M^ r')^(V}(t - r y (r, r'), r')) dr' (1.1) 
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for i = 1, . . . ,p. The intrinsic dynamics exhibit exponential decay with on > for i = 
1, . . . ,p. The propagation delays %(r, r') measure the time it takes for a signal sent by a 
type-j neuron located at position r' to reach a type-i neuron located at position r. For the 
definitions and interpretation of the real valued connectivities Jij and the positive, real 
valued synaptic activation functions Sj appearing in (1.1) we refer to §2 of [50]. 

The aim of this paper is to demonstrate how general theory from the field of delay equa- 
tions can be used successfully to analyse stability and bifurcation of equilibrium solutions 
of (1.1). For this we consider a specific class of delay equations of the form 



where Y is a Banach space, F : C([—h,0];Y) — > Y is a smooth y-valued function on the 
Banach space of continuous y-valued functions equipped with the supremum norm, 



is the history at time t > and <p £ C([—h,0];Y) is an initial condition. The parameter 
h G (0, oo) is a finite delay. As the reader may have expected, the acronym DDE stands 
for delay differential equation. 

Systems of this type naturally extend the case of classical DDE with Y = R n for which 
a rather complete dynamical theory based on perturbative calculus of dual semigroups 
[8], [9], [10], [11], [21] is available in [22]. Recently it was understood that, from an 
abstract viewpoint, various apparently different classes of delay equations can be cast 
and analysed within the same functional analytic framework of dual perturbation theory, 
largely independently of the particulars of a certain class. It is only in the choice of the 
underlying function spaces and the spectral analysis that these details matter. In [18] 
purely functional equations (also called renewal equations) as well as systems of renewal 
equations coupled to delay differential equations are investigated for the M n -valued case 
and finite delay. In [20] the analysis is extended to the case of infinite delay. In [19] 
abstract (Banach space valued) renewal equations with infinite delay are considered. The 
forthcoming paper [47] treats general aspects of abstract equations of the type (DDE). 

The outline of this paper is as follows: in §2 we introduce the functional analytic setting 
and state the equivalence between the abstract delay equation and an abstract integral 
equation using sun-star calculus. We also state a linearization theorem. In §3 we start 
with some general results on the resolvent and spectra, primarily based on [23]. For a 
specific class of connectivity functions, i.e. finite sums of exponentials, and in one spatial 
dimension, we explicitly calculate the spectrum and the resolvent. It turns out that the 
point spectrum is determined by a determinant condition. In §4 we give the normal form 
coefficients for the critical center manifold in case of Hopf and double Hopf bifurcation. 
This is applied in §5 to a scalar neural field equation with a bi-exponential connectivity 
function modelling an inverted Wizard hat. The system is discretised as in [27] and 
the spectrum of the discretised system is compared with the true spectrum, showing 
convergence. We identify in the true spectrum a Hopf point and a double Hopf point. 
For both cases the normal form coefficients are computed, which allows us to identify the 
sub-type of the bifurcation at hand. The theoretical results are confirmed by numerical 
experiments. We end this paper in §6 with conclusions and an outlook on future work. 

Upon finishing this paper we encountered the online preprint [51], addressing similar 
questions. We feel that there are enough substantial differences between the two papers to 
render both of them interesting. Moreover, we have reasons to believe that the choice Y = 
L2(Q) for the spatial state space made in [5 ] leads to non-trivial technical complications, 




(DDE) 



x t (0) := x(t + 6) Vi>0, Oe[-h,0] 
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see §2.4 below. In this paper, we employ sun-star calculus, from which a number of general 
results is immediately available. The center manifold, for instance, was obtained in [22] 
for the abstract integral equation, covering what we need here. 

2. Functional analytic setting 

2.1. Basic definitions and assumptions. It is rather straightforward to associate with 
(1.1) a problem of the type (DDE), but see §2.4. To keep the setting as simple as possible, 
we focus on the single population case p = 1 when (1.1) takes the form 

^(t, r) = -aV(t, r) + J J(r, v')S{V{t - r(r, r'), r')) dv' (2.1) 

For mathematical convenience we extend the spatial domain $7 by its boundary dVt and 
work on Q = S2 U <9Q with Lebesgue measure \Q\ < oo. We formulate a number of basic 
hypotheses on the modelling functions appearing in (2.1). These will be tacitly assumed 
to hold throughout the remainder of this paper. More specific functional forms will be 
chosen in subsequent sections. 

(Hj) The connectivity kernel J G C(S7 x O). 

(Hs) The synaptic activation function S G C°°(M) and its kth derivative is bounded 
for every k G No- 

(H T ) The delay function r G C(£l x $7) is non-negative and not identically zero. 
From (H T ) we see that r is bounded on the compact set 0. Hence we may set 

< h := sup{r(r, r') : r,r'efi}<oo 

Let Y := C(£l) be the Banach space of continuous real- valued functions on with norm 

\\y\\ := sup|y(r)|, 
ren 

We also set X := C([— h, 0]; Y). When (f> G X, t G [— h,0] and r G O we will sometimes 
abuse notation and write <f>(t,r) instead of (fi(t)(r). On X we have the norm 

u\\-.= sup \m, .)|| 

te[~h,o] 

Define the nonlinear operator G : X — > Y by 

G(0)(r) = f J(r,r')S^(-T(r,r'),r'))dr' V^Gl,VrG0 (2.2) 
Jn 

The following lemma is standard, but in light of the difficulties pointed out in §2.4 we 
provide a detailed proof. 

Lemma 1. G : X — > Y is well-defined by (2.2). 
Proof. Obviously, for any (ft € X the map 

[-h, 0] x U 3 (t, r) (ft(t, r) G R (2.3) 

is continuous. 
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Now, given cp € X we consider for points r,r G 17, 

|G(0)(r)-G(0)(r)| < / [J(r, r') - J(r, t>)]S(<K-t(t, r'), r')) dr' 

Jn 

+ f J(r, i>)[SM-t(t, r'),r')) - S(<f>(-r(v, r'), r'))] dr' 
Jn 

<C S I |J(r,r')- J{v,v')\dv' 
Jn 

+Cj f \S(c/>(-r(r, r'), r')) - 5(0(-r(r, r'), r'))| dr' 



where C5 > and Cj > are constants bounding S and J. Let e > be given. By 
the uniform continuity of J on Q x 17 there exists 5/ > such that the first integral does 
not exceed |0|e for all r,f £ $7 satisfying |r — r| < 5j. Regarding the second integral, the 
continuity of (2.3) and (H r ) implies the continuity of 

Tl x II 9 (r, r) ^ 0(-r(r, r), r') € M (2.4) 

Let / C 1 be a compact interval containing the range of (2.4). Then S is uniformly 
continuous on /. Hence there exists 5s > such that \S(u) — S(v)\ < e for all u,v £ I 
satisfying \u — v\ < 5s- Since (2.4) is uniformly continuous, there exists 5' > such that 
|r — f| < 5' implies \<j>(— r(r, r'), r') — 0(— r(r, r'), r')| < 5s for all r' G 17. Consequently, if 
|r — r| < (5 ; then the second integral does not exceed |f2|e. □ 

Using the definition (2.2) of the operator G, we see that studying (2.1) is equivalent to 
analyzing the following initial value problem 

(v(t) = -aV(t) + G(V t ) t>0 
\V(t) = <t>(t) t€[-h,0] 

where V : [—h, 00) — )■ Y is the unknown and <fi G X is the initial condition. Then (NF) is 
of the form (DDE) when we define F : X — > Y by 

F((t>) := -«0(O) + G(<j>) V^GX (2.5) 

with G given by (2.2). The notion of a solution of (DDE), and consequently (NF), is a 
direct generalisation of the solution concept for classical DDE. 

Definition 2. A function x £ C([-h, 00); Y)nC 1 ([0, oo);Y) that satisfies (DDE) is called 
a global solution of (DDE). 

Sometimes we will omit the qualifier global and simply speak of a solution of (DDE). 
We conclude with a simple observation, which follows directly from the fact that (H5) 
implies that S satisfies a global Lipschitz condition. 

Lemma 3. The operator F : X Y defined by (2.5) is globally Lipschitz continuous. 

Proof. It suffices to show that G satisfies a global Lipschitz condition. If tfi, 4> £ X and 
r, r' G 0, then 

|0(-T(r,r'),r') - 4>(-t(t, r'), r')| < sup |<^(-r(r, r), r") - 0(-r(r, r'), r")| 

r"£SJ 

< sup sup \(f)(t,r") -~4>{t,Y") I = 
te[-A,o] r "eii 

Hence we obtain 

||G(0)(r) -G(0)(r)|| < |H| sup J sup S' \\<f> - 4>\\ VrGfi 

nxn R 
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where the suprema are finite due to (Hj) and (Hg). □ 

2.2. Dual semigroups and DDE. In this subsection we provide a very brief introduction 
to sun-star duality and its consequences for the analysis of (NF). For a more complete 
treatment we refer to [22] and, regarding the analysis of abstract DDE, the forthcoming 
paper [47]. 

In this subsection Y will be a Banach space and X := C([— h, 0]; Y). In conjunction 
with (NF) we will assume that Y = C(S7). From an abstract point of view, solving a 
delay equation amounts to obtaining the future state of the system, say at time t > 0, 
from knowledge of the present state. This is done in two steps. First, the present state (a 
continuous function on the time segment [— h, 0]) is extended to the interval [— h, t]. Next 
the part of this extension living on [t — h,t] is shifted back to [— h,0]. Dual perturbation 
theory provides a systematic method to embed X into a bigger Banach space, the so-called 
sun-star dual X *, in which the extension and shifting operations are neatly separated. 
In broad lines, this works as follows. 

If the extension problem is trivial, i.e. if F = in (DDE), then the solution semigroup 
corresponding to (DDE) is the shift semigroup To, defined as 

{T ° m{6) = {Z) d) _ 0<!+r° V*e*.^M€bM (2-6) 

Let Aq be its infinitesimal generator. We represent X* by the space NBV([0, h]; Y*) 
of functions ij : [0, h] —> Y* of bounded variation, normalised such that ry(0) = and 
i](t+) = r](t) for all t G (0,/i). Elements of X and X* are in duality via an abstract 
bilinear Riemann-Stieltjes integral. Since X is not reflexive, the adjoint semigroup Tq 
may not be strongly continuous on X* . Let X® C X* be the maximal subspace of strong 
continuity of Tq . It is easy to see that X is positively Tq -invariant and, moreover, 



X = D(A*) (2.7) 

where Aq is the adjoint of Aq. Let Tq 3 be the strongly continuous semigroup on X 
obtained by restriction of Tq to X . Its infinitesimal generator Aq is precisely the part 
of A* in X , 

D(A°) = {4> & G D(A* ) : A* Q 4> & G X }, A®^ = J$f> 

In [ , Thm. 2.2] it is shown that X may be identified with Y* x T 1 ([0, h];Y*) where 
the second factor is the space of Bochner integrable y*-valued functions on [0,h]. 

Performing this construction once more, but now starting from the strongly continuous 
semigroup T o (t) on the Banach space X , we obtain the adjoint semigroup Tq * on the 
dual space X * and its strongly continuous restriction Tq® to the positively invariant 
subspace X 00 = D(Aq*). The infinitesimal generator of Tq 50 is again given by the 
part of A®* in X 00 . Following [4, §1.2] we suppose that it is not possible to represent 
X * = Y** x [T 1 ([0, h]; Y*)]* in terms of known functions or measures, since Y** does 
not have the Radon-Nikodym property. However, the subspace X 00 of strong continuity 
may be identified with C([— h, 0], Y**), see [29, Thm. 3.11]. Of course this representation 
is only semi-explicit, since a representation for Y** itself is unknown. The original space 
X is canonically embedded into X * via j : X — > X * given by 1 

((/> ,j0) := (<j>,<f> e ) V</> G X, V</> G X (2.8) 



^In this paper we adopt the 'postfix notation' for the action of a functional on a vector. That is, if W 
is a Banach space with dual space W* , w € W and w* £ W* , then (w,w*) := w*(w). 



6 



S.A. VAN GILS, S.G. JANSSENS, YU.A. KUZNETSOV, AND S. VISSER 



Since Y is not reflexive, it follows that the range of j must be a proper subspace of X®®. 
This fact is expressed by saying that X is not sun-reflexive with respect to the shift 
semigroup To, a situation that contrasts the classical case Y = W 1 . 

We proceed to explain how (DDE), and consequently (NF), fits into the above abstract 
context. Define 5 G £(X®, Y*) as 

6<f> Q :=y* V0 = (y*,g) £ X® (2.9) 

Then 5* G £(Y**,X®*). Let £ G £(Y,X®*) be the restriction of 5* to Y, viewed as a 
subspace of Y** . Explicitly, 

(y,6^ G ) = (<p G ,£y) V y G Y, V 0® G X® (2.10) 

Define R : X — > X®* by R := I o F with i 7 as in (DDE). The following lemma will prove 
to be useful in §4.4. Observe that each (y, /) G Y x L°°([— /i, 0]; Y) defines an element of 
X®*. Hence we may identify Y x L°°([—h,0];Y) with a subspace of X®*. 

Lemma 4. R(<fi) = (F(^),0) /or all 4> G X. Hence R maps into Y x {0}. 

Proo/. Let (/> G X and write R(4>) = (y**,w*) G X®* for certain y** G 1"** and w* G 
[LiQC/i]; Y*)]*. Then, for any 0® = (y*,g) G X®, 

(<t>°,R(<l>)) = {y*,y**) + {9,™*) (2-11) 

On the other hand, from the definition of R we obtain 

<0®,ii(0)> = <0®,£F(0)> = (F(0),^®) = <F(0),y*> (2.12) 

where in the second equality we used (2.10) and the third equality is due to (2.9). By 
comparing (2.11) and (2.12) we see that y** acts on y* by point evaluation in F((f>) G Y 
and w* = 0. Hence R(4>) = (F(cj)), 0) and consequently R maps into Y x {0}. 

□ 

Remark 5. In the 'classical' case where Y = W 1 , the previous lemma shows that R is 
a (possibly non-linear) operator of finite rank that takes values in the 'point component' 
Y only, see [22, §§111.3 and VII. 6]. In the present setting with dimY = oo we lose the 
former, but retain the latter property. 

We now consider the so-called abstract integral equation of the form 

u (t) = T (t)cf) + r 1 ^^ T®*(t- s)R(u(s))ds^ Vi>0 (AIE) 

where (p G X is an initial condition, u G C([0, oo); X) is the unknown and the convolution 
integral is of weak* Riemann type, see [22, §111.1] and also [22, Interlude 3.13 in Appendix 
II]. In [47] it is shown that this convolution integral takes values in the range of j. Conse- 
quently, the right-hand side of (AIE) is well-defined. The connection between (DDE) and 
(AIE) is revealed in the following theorem. 

Theorem 6 (Equivalence of (DDE) and (AIE)). Let (f> G X be given and let R = £ o F 
with F G C(X, Y). The following two statements hold. 

(i) Suppose that u G C([0,oo);X) satisfies (AIE). Define x : [— h, oo) — > Y by xq := <fi 
and x(t) = u(t)(0) for t > 0. Then x is a global solution of (DDE) in the sense of 
Definition 2. 

(ii) Conversely, suppose that x is a global solution of (DDE). Define u : [0,oo) — >• X 
by u(t) := xt- Then u G C([0,oo);X) and u satisfies (AIE). 

It is routine to show that (AIE) admits unique global solutions on [0, oo) if F is globally 
Lipschitz continuous. Thus, by Lemma 3 and Theorem 6 we find 
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Corollary 7. For any <fi £ X problem (NF) has a unique global solution. 

Of course, establishing well-posedness for (NF) does not require sun-star duality. Yet, 
it turns out that (AIE) is a very convenient tool in proving many standard results of 
dynamical systems for abstract DDE, such as the principle of linearised (in)stability, the 
existence of stable, unstable, and center manifolds, and theorems on local bifurcation. The 
following linearisation theorem is a direct generalisation of the corresponding result in the 
sun-reflexive case, see [22, §VII.5] or [10]. 

Theorem 8 (Linearisation at a steady state). Let F G C 1 (X, Y) and R = I o F and let 
£ be the strongly continuous non-linear semiflow on X associated with (AIE). Let (f) £ X 
be a steady state o/S, i.e. £(i)(0) = <j) for all t > 0. The following statements are true. 

(i) For each t > the operator E(t) : X — > X is continuously Frechet differentiable in 
(ft with derivative DH(t)((f>) £ C(X). 

(ii) Upon defining T(t) := Z?£ (£)(</>) for each t > one obtains a strongly continuous 
semigroup in C(X). The domain of its generator A is given by 

D(A) = {(p £ X : cp' £ X and cp'(0) = DF(4>)(j)}, Adp = </>' (2.13) 

(hi) For every 4> £ X the function T(-)(j> £ C([0, oo), X) is the unique global solution of 
the linear abstract integral equation 

T{t)<t> = ToWt + j- 1 ^ T *(t - s)£DF($)T(s)<t>ds) 

We observe that the above theorem produces a new strongly continuous semigroup T 
on X with generator A. For this semigroup we may likewise calculate the sun-star duality 
structure, just as we did for the shift semigroup To defined by (2.6). It turns out that the 
spaces X® and, consequently, A * are the same for both semigroups. Indeed, if we put 
B := I o DF{4>) £ C(X, A *) and slightly abuse notation by writing B* £ £(A°, A*) for 
the restriction of the adjoint of B to A®, then just as in the sun-reflexive case [22, §111.2] 
one proves that the adjoint of the generator A of T is given by 

D(A*) = D(A* ), A* = A* Q + B* 

By (2.7) the sun-duals of A with respect to To and T are identical and may both be 
denoted by A®. Moreover, 

D(A e ) = {0° £ D(A*) : A*(f> Q £ A }, A Q = A* (2.14) 

Let A Q * : D(A *) C A * -)• A * be its adjoint. For A G * the situation is slightly more 
difficult than in the sun-reflexive case, because D(A®*) % j(X). (Indeed, if it were true 
that D(A°*) C j(A), then it would follow that A 00 = D(A®*) C j(X) and A would be 
sun-reflexive with respect to To, also see [22, §111.8].) The next lemma is sufficient for our 
purposes in §4.4. 

Lemma 9. If $ £ C 1 ([-/i, 0]; Y) then jcp £ D(A G *) and A G *jcp = (0,0') + (DF($)<f>,0). 

2.3. Differentiability results. The following two results concern the smoothness of the 
operator G defined by (2.2) and appearing in the right-hand side of (NF). When k = 
1,2,... we denote by £fc(A, Y) the space of bounded /c-linear operators from A to Y. 
When k = 1 we write C(X,Y) instead of C\(X, Y). For a review of differentiation in 
Banach spaces, we recommend [3, Ch. 9]. 
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Lemma 10. The operator G : X -4 Y defined by (2.2) is Frechet differ entiable with 
derivative DG(4>) £ C(X, Y) in the point (f> £ X given by 

(DG(M)(r) = [_J(r,r')S'(<P(-T(r,r'),r'm-T(r,r'),r'))dr' (2.15) 
Jn 

for all ip £ X and all r E Q. 

Proof. First we consider the operator DG{4>) denned by the right-hand side of (2.15). 
Using standard methods as in the proof of Lemma 1 one shows that DG((f>)ip £ Y for all 
ip £ X. These steps are omitted. As in the proof of Lemma 3 we begin by noting that if 
r, r' E then 

t(t, r'), r')| < sup \ip(— r(r, r'), r")| < sup sup \ip(t, r")| = \\ip\\ (2-16) 
r"eQ ts[-fe,o] r "en 

This implies that ||(.DGr(0)^)|| < M[|^[| where M > is a constant depending on 0, J 
and 5. Hence DG((j>) £ C(X,Y). 

Next we verify that DG(<p) is indeed the Frechet derivative of G at (p. Introduce the 
shorthand notation (j) T (r, r') := </>(— r(r, r'), r'). For rj £ X and r£S], 

G(0 + 77) (r) - G(0)(r) - [ZX?(<£)r/](r) 

= / J(r, r') [5(^(r, r') + rf (r, r')) - S(0 r (r, r')) - S' (r, r'))^(r, r')] dr' 

Consider the integrand for fixed r'. It follows from the Mean Value Theorem that there 
exists c = c(4>, 77, r, r') E (0, 1) such that 

5(0 T (r,r') +77 r (r,r')) - <S(</> r (r, r')) = ) 7 T (r,r')5'(^(r > r') + erf (r,r')) 

Consequently, 

5(^(r,r') + ^(r,r')) - 5(^(r,r')) - S"(^(r, rO)r? T (r, r') 

= [S'(^(r, r') + of (r, r')) - S'(0 T (r, r'))] if (r, r') 

Since S' is uniformly continuous on compact intervals and \4> T (r, r')| < ||0|| and |rf (r, r')| < 
\\rj\\ for all r, r' E f2, it follows that for every e > there exists <5 > such that 

\S' (<p T (r,r') + crj T (r,r')) - S'((p T (r,r'))\ <e Vr.r'eH 

provided \\rj\\ < 5. Therefore, if 1 1 77 [| < 5 then 

\\G(<l> + T,)-G(tf>)-DG(<f>)T,\\<Me\\ V \\ 

where M > depends on Q and J. This establishes differentiability. □ 

Proposition 11. The operator G defined by (2.2) is in C°°(X, Y). For k = 1,2,... its 

kth Frechet derivative D k G{(p) £ C k (X, Y) in the point <f> £ X is given by 

(Z? fe G(</>)(Vi, • • • , ^fc))(r) = / J(t, r')S^(<f>(-r(r, r'), r')) J] M~r(r, r'),r') dr' 
/or V'l; ■ • • > ipk ^ X and r E 0. 

Proof. For A; = 1 the statement reduces to Lemma 10. Fix k > 2. We need to check 
that D^G : X -4 £fc_i(X,Y) has Frechet derivative D k G{4>) £ C k {X,Y) in the point 
<fi £ X. Again we remark that D k G(<p){ipi, . . . , ip k ) £ Y but we omit the proof. We begin 
by observing that D k G{<p) £ C k (X,Y). Indeed, by (2.16) we have 

\\D k G{4>){i> x ,...^ k )\\ <MU x \\-...-U k \\ 
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for all tpx , . . . , t/jfs; E X, where M > is a constant depending on $7, J and S. 

We conclude by verifying that D k G(<p) is indeed the derivative of D k ~ 1 G at <fi G X. 
Using the same shorthand notation as in the proof of Lemma 10, we consider, for rj £ X, 
ijj = (tpi, • • ■ , ipk-i) e X k ~ x with ||^j|| < 1 for all i = 1, . . . , k — 1 and r G S7, 

(D^G^ + vWir) - (D k - l Gm)(r) - (D k G^)( V ^))(v) 

r fc-i 
= / J(r,r')R(r,r')H^l( ri r') dr ' 

Jn i=1 

where 

R(r,r') := S^i^r, r') + rf (r, r')) - S^if "(r,r')) - ^(^ T (r, r'))?? r (r, r') 

Exactly as in the proof of Lemma 10 we may use the Mean Value Theorem and uniform 
continuity of S^ k ' on compact intervals to conclude that for each e > there exists 5 > 
such that \\rj\\ < 5 implies |i?(r,r')| < e\\r]\\ for all r,r' G O. Hence we have 

WD^G^ + ri)i/> - D k ~ l G{cj))^ - D k G((p)(rj,ip)\\ < Me\\r]\\ 

provided ||t/|| < 5, where M > depends on £1 and J. □ 

2.4. Choosing the spatial state space. We believe that our choice for Y = C(fi) 
made in §2.1 deserves some comments. In [27, 50] the authors instead elect to work 
with the Hilbert space Y = L 2 (Q). In our opinion this choice suffers from at least three 
mathematical complications. 

The definition of G. It is no longer clear that G is well-defined by (2.2). Namely, apart 
from square integrability one also needs to verify the following. If (f>, <f> £ X and for all 
t G [—h, 0] one has 

(j>(t,r') =4>(t,r') a.e. r' G Q (2.17) 

(where a.e. stand for almost everywhere, i.e. 4>{t, •) and 4>(t, •) represent the same element 
in L 2 (Q)) then this should imply that for almost all r G Q one has 

0(-r(r, r), r) = 0(-r(r, r), r') a.e. r G ft (2.18) 

There are bounded r G C(J7 x $7) for which this implication fails. For example, let 
(7 = (0, 1), write x = r and r = r' and let -0 and ^ be representatives of the same element 
m L 2 (R) that differ in zero. If we define 

<j){t,r) := Tp(r + t), 4>(t, r) := ^(r + t) Vr G S7 

then (p,(j) £ X and (2.17) holds for all t G [—1,0]. However, if r(x,r) = r, independent 
of x G S7, then 0(— r(x, r), r) = ^(0) and 0(— r(x, r), r) = ^(0) which shows that for all 
x = r G O we have inequality in (2.18). Clearly, this choice of r is very unrealistic, but 
it does indicate a problem that needs to be addressed when one works with spaces of 
equivalence classes of measurable functions. It is not obvious that a more realistic choice 
such as r(x,r) := \x — r\ does not exhibit the above phenomenon. 

First order Frechet differentiability. Even if we assume that the above problem can be 
solved satisfactorily by imposing additional (physiologically plausible) conditions on r, 
there remains the question of whether the first order Frechet derivative of G appearing in 
Lemma 10 maps X into Y when Y = L 2 (f2). For the sake of simplicity, let us assume that 
J(r, r') = 1 and = 0. Then the mapping 

Sl3r4 f ■;/;(— r(r, r'), r') dr' G K (2.19) 
Jn 
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should be in L 2 (0) for all ift G X. This is not obvious. An attempt to prove this statement 
is contained in the proof following [_•. r, Lemma 3.1.1]. The authors write, for r G fi, 



The first estimate is by the Cauchy-Schwarz inequality. As it stands, the second estimate 
only seems to be valid under certain extra conditions on ift and / or r, since r depends on 
the integration variable. 

Higher order Frechet differentiability. In verifying second order differentiability we en- 
counter problems similar to those pointed out above. Another complication appears in 
conjunction with derivatives of order three and higher. For instance, consider k = 3 
in Proposition 11. Let 0, = (0, 1) and write x and y for r and r'. Define if G X by 
ift(t, r) := r~5 for all t G [—h, 0] and r G $7. Then clearly if) G X but the integral 



Which space to choose instead? It appears that the choice Y = L P (Q) with 1 < p < 1 is 
not very fortunate. Moreover, from a biological point of view it is rather unclear why the 
membrane potentials should be merely p-integrable on and not necessarily bounded. 

Thus we are led to consider alternatives. Within the class of Hilbert spaces the Sobolev 
space H k (Q.) comes to mind. By standard Sobolev embedding theory each element of 
H h (Q) has a (unique) continuous representative, provided k G N is sufficiently large 
(depending on the dimension of Q). Moreover, H k (Q) is a Banach algebra under mild 
conditions on Q [1, Thm. 5.23]. However, for arbitrary (ft G X the mapping (2.19) cannot 
be expected to possess k weak derivatives in L 2 (J7). 

Other possibilities are Y = L°°(fi), Y = B(H) and Y = C(H), where B(H) is the 
Banach space of everywhere bounded, measurable functions on Q. Note that the first two 
spaces differ in the sense that L°°($7) consists of equivalence classes of essentially bounded, 
measurable functions on f2. The first choice satisfies all our needs, but it may potentially 
suffer from the problem indicated in §2.4. The second choice takes care of all the above 
technical complications but also introduces new ones. Most notably, the Arzela-Ascoli 
theorem, used in §3.1, does not hold in B(Q). The choice Y = C(O) seems to be fitting 
both from a modelling as well as from a technical perspective. 





diverges for all x G so D 3 G(0) does not map X into Y. 



3. Resolvents and spectra 



Let DG{4>) G C(X, Y) be the Frechet derivative of G at the steady state vector (ft G X, 
i.e. (ft is independent of time (but possibly dependent on space) and 

-a(ft + G(j>) = (3.1) 
by (NF). Using Lemma 10 we obtain 





(3.2) 



where 



J (r,r') := J(r,r')S'(0(-r(r,r'),rO) 



(3.3) 
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In this section we are interested in the spectral properties of the linear problem 

V(t) = -aV(t) + DG(4>)V t t > 
V(t) = <f>(t) t€[-h,0] 

with a > 0, which is a special case of the problem 

x(t) = -ax(t) + Lx t t>0 
x(t) = cj){t) t G [-h,0] 

where Y is a complex Banach space, X = C([— h, 0]; Y), L G C(X,Y) and a G C 

Remark 12. For the spectral analysis of this section it is necessary to work in Banach 
spaces over C. So, whenever we discuss the spectral properties of (3.4), we implicitly 
assume that the spaces X and Y and the operators acting between them have been com- 
plexified. In fact, one should also complexify the sun-star duality structure introduced in 
§2.2. This task is not entirely trivial and rather tedious. Fortunately it has been carried 
out in [2! , §111.7]. 

In §3.1 we make several standard observations on the structure of the spectrum of the 
generator of the strongly continuous semigroup solving (3.5). Familiarity with the basics 
of spectral theory and semigroup theory is presumed, for which we recommend [45, Ch.V] 
and [23]. We would also like to mention the nice application- inspired paper [5] for a 
detailed treatment of abstract linear DDE with bounded right-hand sides, partially in the 
context of Hale's [ ' ] formal duality approach. Some of our statements are similar to those 
found in [50, §3.1], but our approach (as well as the choice of state space, see the remarks 
in §2.4) is sometimes different. For instance, following [22, Def. II. 4. 22] and [ , §4.1] we 
believe that the employment of Browder's (instead of Kato's) definition of the essential 
spectrum leads to somewhat simpler arguments. 

In §3.2 we specialise to (3.4) and choose Y = C{Q) and L = DG{<p). It is shown how 
to obtain explicit representations of resolvents and eigenvectors for a particular (but still 
rather general) choice of connectivity function J. 

3.1. Spectral structure. We recall from Theorem 8 in §2.2 that the strongly continuous 
semigroup T on X corresponding to the global solution of (3.5) is generated by A : D(A) C 
X — > X where 

D(A) = {(f) G X : (j)' G X and 0'(O) = -a<f)(0) + L<f>}, A0 = <\> (3.6) 

At this point we establish some standard notation. Let S : D(S) C U — > U be a closed 
linear operator on a complex Banach space U. We denote by p(S) C C, cr(A) and cr p (A) 
the resolvent set, the spectrum and the point spectrum of S, respectively. When z G p(S) 
we write 2 R(z, S) := [z — S)^ 1 for the resolvent of S at z. For any z G C we let 1Z(z — S) 
and M{z — S) denote the range and the nullspace of z — S. 

The results in this subsection are rather easy consequences of the following generalisation 
of [22, Thm. IV. 3.1 and Cor. IV. 3. 3]. It will turn out to be very convenient to employ 
tensor product <S> notation as introduced in [23, p. 520]. We recall the definition from 
there for the reader's convenience. 

Definition 13. Let U, V be complex Banach spaces and let J- (I, V) be a complex Banach 
space of 1^-valued functions defined on an interval ICR, Let B G C(U, V) and g : I —> C. 



It is customary to suppress the identity operator and write A — S instead of XI — S. 



(3.4) 
(3.5) 
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If the map g (g) v : I 3 s \— > g(s)v G V is in J- (I, V) for all v G V, then we define 
g®B:U -»■ ^(J, F) by 

[(5 ® B)«](a) := (<7 ® Bu)(s) = g(s)Bu 
for all it G C7 and s £ I. 

We also introduce some auxiliary operators. For each z G C and # G [— h, 0] we set 
£z(0) : = e ze - With L as in (3.5) we define 

L z eL{Y), L z f:=L(e z ®f) (3.7) 

^ G C(X), (H Z <P)(9) := / e^- s U(s)ds 

Je 

S z eC(X,Y), S z (f> := (t>(0) + LH z <p 

for all / G Y, (p G X and G 0]. 

Proposition 14 ([23, Prop. VI. 6. 7]). For every z G C define A(z) G £(y) 6y 

A(z):=z + q-L 2 (3.8) 

Tften (j) G 72- (2; — A) if and only if 

A(z)f = S z <f> (3.9) 

has a solution f G Y and, moreover, z G p(A) if and only if f is also unique. If such is 
the case, then 

R(z, A)<f> = {e z &{z)- l )S z (t) + H z cj) (3.10) 

Furthermore, S z is surjective for every z G C, so A G o~(A) if and only if G <r(A(A)). 
Finally, tp G D(A) is an eigenvector corresponding to A if and only if ip = E\ ® q where 
q G Y satisfies A(X)q = 0. 

Corollary 15. Let z 7^ —a. If L z is compact, then TZ(z — A) is closed. 

Proof. From the part of Proposition 14 regarding (3.9) we have eft G 1Z(z — A) if and 
only if S z (j) G TZ(A(z)). From the theory of compact operators [45, §5.5] it follows that 
A(z) = z + a — L z has closed range, since z + a ^ 0. Now let (4> n )n&N be a sequence in 
1Z(z — A) converging to some <p G X. Then the sequence (S , z </> n )neN in 7£(A(z)) converges 
to S z cj) G K(A(z)), since TZ(A(z)) is closed. Hence <f> G TZ(z - A). □ 

Remark 16. For the particular case (3.4) with Y = C(Q) and L = DG(cp), compactness 
of L z for each z£C follows easily from the Arzela-Ascoli theorem since L z is a Fredholm 
integral operator with continuous kernel J§e~ ZT . 

As dimY = 00 the shift semigroup Tq on X is no longer eventually compact. Conse- 
quently we need to consider the possibility that a(A) contains points that are not isolated 
eigenvalues of finite type. 

Definition 17 ([7, Def. 11]). The Browder essential spectrum o" ess (S') of a closed and 
densely defined operator 5 : D(S) C U — > XJ consists of all A G a(S) for which at least 
one of the following three conditions holds: 

(i) A is an accumulation point of cr(S); 

(ii) TZ(X — S) is not closed; 

(iii) Ufc>o-^[(^ — S) k ] has infinite dimension. 



ON LOCAL BIFURCATIONS IN NEURAL FIELD MODELS WITH TRANSMISSION DELAYS 13 



We also recall that if A is in the point spectrum cr p (S), then the closure of the subspace 
appearing in (iii) is the generalised eigenspace M(A, S) corresponding to A. Its dimension 
m\ (which may be oo) is the algebraic multiplicity of A. If m\ < oo then A is called an 
eigenvalue of finite type. If m\ = 1 then A is called a simple eigenvalue. 

Corollary 18. Suppose L z is compact for all z ^ —a. Then a ess (A) C {—a}. Moreover, 
cr (A) \ {—a} consists of poles of R(-,A). The order of X as a pole of R(-,A) equals the 
order of zero as a pole of R(-, A(A)). 

Proof. Let A G o~(A) and A ^ —a. Corollary 15 implies that (ii) in Definition 17 (with 
S = A and U = X) cannot be true. Since A + a ^ is in cr(L\) and L\ is compact, 
it follows (again from general spectral theory, see e.g. [45, §5.8]) that A + a is a pole of 
R(-,L\), say of order k > 1. If we can prove that A is a pole of order k of R(-,A), then 
we are done. Indeed, it then follows that in particular A is isolated in o~(A), so (i) in 
Definition 17 cannot hold. By [45, Thm.5.8-A] the same is true for (iii). 

Let us therefore prove that A is a pole of order k of R(-, A). First we remark that the 
map 

C 3 z i-> L z G C(Y) (3.11) 

is continuous at A. The proof of this fact is standard and has been omitted. If z is in p(A) 
then z + a G p(L z ) and 

A(z)- 1 = (z + a- L z y l =[ z + a -{L x + {L z - La))]" 1 

A continuity property of the resolvent [36, Theorem IV. 3. 15] together with the continuity 
of (3.11) at A implies that for z sufficiently close to A we have z + a £ p{L\) and 

A(z) _1 = (z + a- L x )~ l + o(\X- z\) as z -4 A 

where o(|A — z\) denotes a term that vanishes as z — > A. By (3.10) we see that for z 
sufficiently close to A, 

R(z,A) = {e z (z + a - L X )- 1 )S Z + H z + o{\\ - z\) (3.12) 

where it was used that \\S Z \\ remains bounded as z — > A, which can easily been seen from 
(3.7) (with A replaced by z). This already establishes that A is an isolated singularity of 
R(-,A). To conclude the proof we recall that A + a is a pole of R(-,L\) of order k > 1. 
Hence A itself is a pole of order k of the mapping 

C3 z^ (z + a- Lx)- 1 €C(Y) 

The result now follows from (3.12) since C3z4 H z £ C{X) is analytic in z = A and the 
zero-order term in the power series expansion of C 3 z i->- S z G C(X, Y) at z = A does not 
vanish, as is easily checked. 

□ 

Hence, although for the application to (3.4) that we have in mind essential spectrum 
exists in the form of the exceptional point —a < 0, it is properly contained in the left half- 
plane and therefore rather harmless. This situation seems to be quite common in DDE 
arising in population dynamics, see the remark in [5, p. 321]. As a pleasant consequence, 
most of the results in [22, §IV.2] have immediate analogues in the present setting. We 
will limit ourselves to the statement of two such results that are also important for the 
application of center manifold theory in §4. 

Lemma 19 ([23, Thm. VI. 6. 6 and Corollary IV.3.11]). The semigroup T generated by A 
is norm continuous for t > h. Consequently ujq = s(A), where ojq is the growth bound of 
T and s{A) is the spectral bound of A. 
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The above lemma implies that, for the linear problem (3.4), the (in)stability of the 
zero solution may be inferred from the location of the poles of R(-,A) in the complex 
plane. More precisely, we have the following result, which is a direct analogue of [22, 
Thm. IV.2.9]. 

Proposition 20. Suppose (3 > —a. Let 

A = A(/3) := {A G a(A) : ReX > (3} 

and let P A G C(X) be the spectral projection associated with A, see [45, §5.7]. Then 

X = 1Z(P A ) e K{I - P A ) 

where the first summand is finite dimensional, the second summand is closed and both 
summands are positively T -invariant. Moreover, there exist K > and e > such that 

\\T(t)P A \\ < Ke^ +£)t \\P A \\ Vt < (3.13) 

\\T(t)(I - P A )\\ < Ke (f3+£)t \\I - P A \\ Vt > 

We observe that T(t)P A is well-defined in (3.13) for all t < 0, since T(t) extends uniquely 
to a group on the finite-dimensional range of Pa- 

The extension of the above decomposition and exponential estimates to X * proceeds 
exactly as in [ , p. 100 - 101]. 

3.2. Explicit computations. In the remainder of this section we consider a homogeneous 
neural field with transmission delays due to a finite propagation speed of action potentials 
as well as a finite, fixed delay To > caused by synaptic processes. Space and time are 
each rescaled such that ft = [— 1, 1] and the propagation speed is 1. This yields 

t(x, r) = To + \x — r\ Vx,rGf2 (3-14) 

For the connectivity function we take a linear combination of N > 1 exponentials, 

N 

J(x,r) = J2^e~^ lx ~ r{ Vi.reH (3.15) 
i=i 

where 

Ci G C with Cj ^ 0, fii G C with yUj / fij for i / j 

(As the number N of exponentials remains fixed, we suppress it in our notation.) In 
addition to (H5) we also require here that S(0) = 0. We study the stability of a spatially 
homogeneous steady state eft = by analysing the spectrum of the linearised system (3.4). 
Following (3.3) we incorporate S'(0) into the connectivity function, 

N 

J (x,r) = ^ Ci e-^- r l, c i = S'(0)c i 
i=i 

In order to avoid overly convoluted notation we henceforth write J instead of Jo- Assuming 
the form (3.15), in the next two subsections we explicitly compute the point spectrum 
o- p (A) with A as in (3.6) as well as the resolvent operator P(A, A) for A G p(A). 
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3.3. Characteristic equation. In this example the operator A(A) introduced in (3.8) is 
given by 

(A(X)q)(x) = (X + a)q(x) - J J(x, r)e~ XT ° e~ X \ x - r \q{r) dr (3.16) 

for all A G C, q G Y and x G U = [-1, 1]. We let 

ki-.= X + tM Vi = l,...,JV (3.17) 
and define for each i = 1, . . . , N the integral operator Ki G CiY) by 

(K ig )(x) = j e- k ' lx - rl q{r)dr 

and set (Kq)(x) : = [(K\q)(x), . . . , (KNq)(x)] G C^. By introducing c := [ci,...,cjv] G 
C , A (A) is written as 

A(X)q= (X + a)e Xr °q- (c-Kq), (c • Kq)(x) := (c • Kq(x)) (3.18) 

where (a ■ b) := YliLi a i^i 1S a pairing of two complex vectors a = [ai, . . . , a n ] and 6 = 
[bi, . . . ,b n ]. We solve the equation A(X)q = by formulating a linear ODE in terms of q 
by repetitive differentiation. For this purpose the next lemma is useful. 

Proposition 21. All solutions q £Y of the equation A(X)q = are in fact in C°°(0,). 

Proof The range of Ki is contained in C 1 (O) and therefore any solution of the equation 
A(X)q = is an element of C°°(fT). □ 

Let q G C 2 (0). The first derivative of A(X)q with respect to the spatial variable contains 
terms that involve integration over the intervals [— l,x] and [x, 1]. The second derivative 
has a nicer structure: 

D 2 x A(X)q = (A + a)e XT0 q (2) + 2(c • k)q - (ck 2 ■ Kq) (3.19) 

in which q^ denotes the second derivative of q and for each m G N the vectors k m and 
ck m in have elements fef 1 and Cikf 1 , respectively, for i = 1, . . . , N. This identity allows 
for straightforward calculation of higher derivatives. For the following lemma we recall 
the definition in (3.17). 

Lemma 22. The set S := {A G C : 3i,j G {1, . . . , N},i ^ j, s.t. k 2 = k 2 } contains at 
most t}N(N — 1) elements. 

Proof All ki are distinct since (by definition) all [ii are distinct. So for % ^ j, k 2 = k 2 
A = — + /ij) G 5. The number of (unique) elements in this set is at most the number 
of unique pairs (i,j), i, j < N,i ^ j, which equals \N{N — 1). □ 

Lemma 23. Let X £ S. Then there exist unique vectors £ = [Co, • • • , Gv-l] £ & N an d 
(3 = [/3q, . . . , /3n] G C n+1 , depending on X and such that for every q G C 2N (Q) one has 

(Co + CiD 2 x + ... + (n-iD 2 x n - 2 + D 2 x N )A(X)q = ((3 + $\D\ + ... + (3n-iD 2 x n ~ 2 + (3 N D 2 x N )q 

Proof. Let Q := [q,q^ 2 \ . . ■ ,q^ 2N ^]. Repeated differentiation of (3.19) yields the following 
system of equations: 

A(X)q 
D 2 x A(X)q 



D™A(X)q 



MQ - V (3.20) 



1(5 



S.A. VAN GILS, S.G. JANSSENS, YU.A. KUZNETSOV, AND S. VISSER 



where 











... 0" 




(c 


Kq) 




(c-k) 





... 




(ck 2 


Kq) 


M := e Ar °(A+a)/+2 


(c-k 3 ) (c-k) 





... 




(ck 4 


Kq) 




(c-k™- 1 ) ... 


(c • A; 3 ) 


(c-k) 




_ (ck 2N 


Kq) 



and I is the identity matrix of size N + 1. (Note that (3.20) is an equality that holds 
on fi. Also, the definition of H T is not used in the current proof, but will reoccur in 
Appendix A.) We take a linear combination of the rows in (3.20) with the components of 
the vector Z := [£, 1] G C^ +1 such that this combination of the elements of V in (3.20) 
vanishes. This eliminates all integral terms entering (3.20) via Kq. Thus we seek £ such 
that Z T V = 0, i.e. 



[C i] 



1 

kl 
ki 



k\ 



2N 



k 2 
k 2 



k 2N 
fi, 2 



k 2 
k 4 



k 2N 
N 



c\K\q 
C2K 2 q 

c N K N q 



(3.21) 



-~w T 



on J7. If this equation is to be satisfied for any q, then we must have WZ 
equivalent to 



1 



k 2 

N 



1 
^2 



k A 



2N-2- 



k 



2N-2 



■ 2N-2 



k\ N 
k 2N 



k 2N 



0, which is 



(3.22) 



it 



The N x N Vandermonde matrix W is invertible since all kf are distinct by Proposition 
22. Hence ( can be found by applying W~ x to (3.22). To find we apply the row vector 
[C, 1] T from the left to (3.20) to infer that /3 T = [£, l] T M. Hence 

~kl N ~ 



(3 = M 



-M 1 



w- 1 






1 



k 2N 



k 2N 

K N 



(3.23) 



□ 



which concludes the proof. 

Remark 24. For A € S the vectors £ and /3 still exist, but they are not unique, as can 
be seen from (3.22). For simplicity we do not consider this case here. 

Theorem 25. Suppose A S and let {(3i}f =1 as in Lemma 23. Then A(A)g = implies 
fog + Piq {2) + ■-- + (3 N -iq {2N - 2) + PnQ^ 







(3.24) 



Proof. Since A(X)q = on O it holds that D m A(X)q 
follows from Lemma 23. 



for all m G N. The result now 

□ 
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Our next objective is to obtain what one could call a converse to the above theorem. 
Specifically, we ask when for a given A G C with A G" S a solution q of (3.24) also satisfies 
A(A)g = 0. For this we start by noting that eigenvalues of the ODE (3.24) are roots of 
the characteristic polynomial 



V(p) = I3 NP 2N + ^jv-ip 2JV - 2 + • • . + fa 1 + A) 



27V-2 



(3.25) 



Evaluating the coefficients (3i of this polynomial by means of (3.23) yields the following 
result. Its proof may be found in Appendix A. 



Proposition 26. For A S the characteristic polynomial V is given by 

Ar M , \ N N N 

v{ P ) = e (A 9 +a) U(p 2 - w) 2 ) + E «hw n ^ 2 - k ^ 2 



(3.26) 



i=l 



i=i 



Since "P is an even function, it follows that if p G C is an eigenvalue of (3.24) then the 
same is true for —p. 

Proposition 27. If (3.26) has 2N distinct roots ±pi(A), . . . , ±/?at(A) then the general 
solution of (3.24) is of the form 



N 



(s) = £[7 i e«W* + 7_ ie - 



(3.27) 



where the coefficients 7-y G C are arbitrary. 

For (3.27) to satisfy A(A)(/a = 0, from (3.18) we see that 



A? 



= (A(A)g)(x) = e Ar °(A + a) ]T [7^ + j-i.e^] 



7=1 i=l 



(3.28) 



e -fcj|x-r|+p<r ^ r _|_ 



7-i 



-fcj|a;-r|-pir ^ r 



must hold for all a; G f2. For notational convenience we have suppressed the dependence 
on A of q, p% and k. Recalling that = [— l,x] U [x, 1] for each fixed 1 £ O, we split the 
domains of integration accordingly. If 



k j (X)^± Pl (X) = 1,2,..., N 



(3.29) 



then (3.28) becomes 

N 

= e Aro (A + a) [li^ X + 7-*e 

i 

N N 



h 2 - n 2 
j=l i=l ^3 p i 



qPi X Tl ' {Mj x 



kj pi 



e -{kj+pi) 

1 g k j x 

kj + Pi 



N N 

3=1 i=i 



2 k 



k 2 - n 2 
K i Pi 



'0 e ~Pi x 



e -{kj+pi) e ~(kj-pi) 



kj + pi 



kj pi 



e J 
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Sorting the terms according to their exponents in x while again suppressing dependence 
on A of pi and k yields 







N N 

]T{7^[e Aro (A + a )-]r 



i=\ 
N 



=1 K 3 Pi 



N N 



' kj - pi 



£ 



7-) 



+ 7-i e " 



kj + pi 



A 2c k- n 



N 



i=l ^ 



£ 



7-i 



Proposition 26 guarantees that all coefficients of e^P*w x vanish. As for the remaining 
terms, all coefficients of e ±fe J^ a! should vanish as well. Thus we must have 



N 



N 



£ c ; e ~ fej ' e ^[E^^— +E 

i=l i=l J ri i=l 

N N 



7-i 



fej + Pi 



E c i e " %e ~^[E^X^£: + E^ 

1=1 



' kj + pi 



i=l 



kj Pi 



where dependence on A of pi and k was suppressed. This yields a set of 2N linear equations: 
one for each e ± k iW x m With T = [71,72, • • • , 7tv , 7-1, 7-2, • • • , 7-jv] and the matrix 5(A) 
defined by 



5(A) := 



(3.30) 



where 



[5 



=Pi(A) 



A JJ,» 



A + /ij -pi(A) ; 



[5 



-«(A) 



A JJ,« 



A + fij + pi(X) 



we seek T such that 



5(A)r = 



(3.31) 



In order for this system to have a non-trivial solution V = T\, it is necessary (and sufficient) 
for the determinant of S(X) to vanish, 



det 5(A) = 

This result is summarised in the following theorem. 



(3.32) 



Theorem 28. Suppose that A 5 and assume that the characteristic polynomial V in 

(3.26) has 2N distinct roots, denoted by ±pi(A) for i = 1, 2, . . . , N. If X satisfies (3.32) 
and (3.29) then A E a p (A). The corresponding eigenfunction is E\ (8) q\, with q\ given by 

(3.27) with T\ a solution of (3.31). 

Remark 29. Two comments on the above Theorem seem in order. 

(i) The above procedure can easily be adapted to cover the degenerate cases excluded 
in Theorem 28. All we need is to adjust the form of q\ in Proposition 27. We do 
not pursue this for reasons of clarity and readability. Rather, in specific instances we 
check that degeneracy is not an issue. 

(ii) We expect that that the order of A as a root of (3.32) equals the multiplicity of A as 
a pole of R(-, A), see Corollary 18 in §3.1. This would give an explicit way to verify 
simplicity of critical eigenvalues in §4. We intend to comment on this issue in future 
work. 



ON LOCAL BIFURCATIONS IN NEURAL FIELD MODELS WITH TRANSMISSION DELAYS 19 



3.4. Resolvent. Now that we are able to reduce determining the point spectrum, in 
this specific example and modulo a technical restriction, to a finite dimensional matrix 
problem, the next step is to determine the solution of the resolvent problem, 

(z - A)ip = <p (3.33) 

i.e. to find a representation of E X in terms of the given function if el when z S p(A). 
For this task we see from Proposition 14 in §3.1 that we first need to solve 

A(z)V(0) = S z <t> (3.34) 

For our specific example the above is equivalent to an integral equation for q := ip(0), 



(z + a)q(x)- I J(x,r)e- ZT0 ~ zlx - rl q(r)dr = h z (x) ViGfi 



(3.35) 



where 



h g (x) := 0(0, x) + 



1 r0 



J(x,r)e- z( - T0+s) - zlx - rl (t>(s,r)dsdr (3.36) 

-1 J— to— \x— r\ 

for all x £ S7. Inspired by (3.27) we propose the following variation-of-constants Ansatz 
for its solution 

N 

q(x) = g(x) + ^2 [li{x)e piX + l-i{x)e~ PlX ] Vx € U 

i=l 

where p±i{z) are distinct roots of (3.25). We seek g G C(O) and 7-1-1, . . . , 7±jv G C 1 (il). 
Substitution into (3.35) and suppressing dependence on z of h, pi and k yields 

v 



>h(x) =e ZT °{z + a)g(x) + e ZT °(z + a) ^ [^{x)e p ' x + 7 _ I (a;)e- p ' 2: ] 



JV 

E- 

3=1 

N 



i=l 
N 



e- k i r g(r) dr + E / Hi r )^ k]+Pl)r + J-i(r)e^- pi >] rfr} 



- Cje.~ kiX { / ^ r g{r) dr + ^ / [ ll {r)e^ +p ^ r + ^ t {r)e (k ^ p ^ r ] dr\ 

j=l i=l 

If (3.29) holds, we may integrate by parts and rearrange the terms, 

JV 

e ZT °h(x) =e ZT0 (z + a)g(x) + e ZT °{z + a) ^ [^{x)e piX + -f_i(x)e- piX ] 



1=1 



-PiX 



o PiX 



- PiX 



N 



JV 



JV JV 



e -fei+Pi 7 .(l) ^ e - fe J-«7_i(l) 
kj Pi kj + pi 

JV 



-kjr 



v-^ e Pir e~ p * r ] 1 



JV JV 

E^1E[- 



-fcj-P* 7i (_l) e -fc*+w 7 _i(_i) 



fcj + Pi 



jv 



E e P,r e -~Pir 
77^ ^'( r ) + I 7 - 



-i(r)]dr} 



(3.37) 
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where again dependency of h, pi and k on z was suppressed. When z G" S, Proposition 26 
is applied and all terms involving e ±Pl ^ x drop out. We can choose g = g z as 

/ \ h z{x) 

9z{x) 



z + a 



provided we can achieve that the remaining terms (i.e. the last four lines) of (3.37) vanish. 
So for j = 1, 2, . . . , N it should hold that for every x G Q, 



'"{.9M + E[ 



r* , N r 

y i ^ r {- 3 (r) + ^[ 



fej — p 



fc 3 + 



— 7-iMl } dr - e fej E [■ 



fci + Pi 



fcj - p. 



y^(r)]}dr + e- fe ^[ 



= o 



(3.38) 



with the same notational convention as before. We seek functions j±i such that the 
integrands and the remaining terms in (3.38) vanish. This yields the system 



^ z M z 

z 



P+(x) 
P-{x) 



T'{x) 



h z (x) 
z + a 



Vx G Q 



where 1 G M. is the vector with one on each entry, 

1 



kj(z) ± pi{z)' 



Pf{x) :=^ N [^'>,...,e 



±p N (z)xl 



(3.39) 



and 



T = [71, • • . ,7jv,7-i, • • -,7-iv] 
If the matrix T(z) is invertible, we find r = T z by matrix inversion and integration, 



r 2 (x) = r , 2 + 



x h z (r) 



X 



z + a 



P~{r) 
P+(r) 



T(z. 



-1 



dr 



(3.40) 



:=T z (x) 

for some initial reference point xq in O and integration constants ro )2 G C 2Ar . Any choice 
of integration constants results in a choince for T z for which the integral terms in (3.38) 
vanish. In order to satisfy the remaining terms in (3.38), Tq <z is chosen as 



0,2 



-S(z)- 1 



s~ sf 

S+ S7 



r,(i) 

TJ-l) 



(3.41) 



for S(z), St, and S 1 " as in (3.30). Clearly, S(z)~ 1 exists if and only if det S(z) ^ 0, 
which is consistent with the fact that the resolvent operator R(z, A) is not defined when 
z G (Jp(A). We are now ready to formulate the key result of this section. 

Theorem 30. Suppose that z G p(A) and 

• z G" S; 

• the characteristic polynomial V has 2N distinct roots; 

• condition (3.29) holds and; 

• the matrix T{z) is invertible. 

Then the solution of (3.33) is given by ip z = e z % q z + H z (j) with 

N 



z + a 



Vx G O 



(3.42) 



i=l 



u>ii/i T 2 given by (3.40) and h z is as in (3.36). 
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Remark 31. The fourth condition in the above Theorem seems peculiar and of a different 
nature than the first three, which already occured as simplifying conditions in §3.3. We 
refrain from investigating this issue here. Wherever we need the result of this theorem, 
we check the fourth condition explicitly. 

4. Normal forms for local bifurcations 

Let 4> £ X be a stationary point of the semiflow generated by (DDE). By Theorem 8 
in §2.2 the linearisation of this semiflow at <p defines a strongly continuous semigroup T 
of bounded linear operators on X , generated by A as in (2.13) and (3.6). If F is as in 
(2.5) then T will be the solution semigroup of the linear problem (3.4), which is of the 
form (3.5). In the present section we prepare for the computation of a critical normal 
form when (j> undergoes a Hopf or a double Hopf bifurcation. The actual computation is 
performed in §5. 

In §2.2 we alluded to the fact that a reformulation of equations of type (DDE), such as 
(NF), as an abstract integral equation of type (AIE) allows for a relatively straightforward 
application of basic dynamical results such as the center manifold theorem. Indeed, by 
(AIE) and the exponential estimates of Proposition 20 in §3.1 the general center manifold 
theory for AIE presented in [22, Ch. IX] is directly applicable to (DDE) in the setting of 
§2.2. We shall relegate a more detailed technical presentation to the forthcoming paper 

There exists an efficient approach based on Fredholm's solvability condition towards the 
derivation of explicit formulas for critical normal form coefficients of local bifurcations of 
dynamical systems. Once such formulas have been derived for a certain class of dynamical 
systems, they may be evaluated for specific equations using spectral information from 
the linearisation at the critical equilibrium or fixed point, together with information on 
the higher order derivatives of the particular non-linearity. The technique goes back to 
[16] and has been successfully applied to ordinary differential equations [37], [38, §8.7]) 
and iterated maps [42], [40], [28]. The resulting formulas have been implemented in the 
software packages CONTENT [39], its successor MATCONT [ ] and CL_MATC0NT for maps. 

In the forthcoming paper [35] the method is applied to AIE and DDE. Here we briefly 
summarise the results related to Hopf and double Hopf bifurcations, obtained using the 
Fredholm solvability technique, see Lemma 33 below. In the Hopf case the corresponding 
formulae have been first obtained in [46, 22] using a different method. As expected, the 
formulae given below look very similar to those given in [37] and [38, §8.7]. However, 
one should pay special attention to their proper interpretation in the current functional 
analytic context. 

4.1. Preliminaries. In §§4.2 and 4.3 we will consider the situation that (j) £ X is a 
stationary point of the non-linear semiflow generated by (DDE) and the linearised problem 
takes the form (3.5) with A as in (3.6) and L z compact for all z ^ —a. There is no loss of 
generality in assuming that = 0. Suppose that A has n c > 1 simple eigenvalues on the 
imaginary axis, counting multiplicities. 

Remark 32. One may show that a (A) = a (A*) = a(A G ) = o-(^ *), see [22, p. 100 - 101] 
and also [23, Proposition IV. 2. 18]. We will use this fact in the remainder of this section. 
For a detailed discussion of the 'lifting' of the spectral properties of A to corresponding 
properties of the various (adjoint) generators, we refer to [22, p. 100 - 101]. 

This implies the existence of a non-trivial center subspace Xq of finite dimension n c 
and spanned by some basis <3? consisting of (generalized) eigenvectors corresponding to 
the critical eigenvalues of A. There exists a locally invariant center manifold Wf oc that is 
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tangent to Xo at the origin. One can show that on Wrl the solution satisfies the abstract 
ODE 

u(t) = j' 1 (A G *ju(t) + R(u(t))) VtSl 

where the non-linearity R is given by Lemma 4 and is as smooth as the mapping F 
appearing in (2.5). Let be the projection of u(t) onto Xo. Then can be expressed 
uniquely relatively to The corresponding coordinate vector z{t) of satisfies some 
ODE that is smoothly equivalent to the normal form 

z{t) = 9vz v (t) + 0(\z{t)\ i ) VtGR (4.1) 

1<M<3 

with unknown critical normal form coefficients g v £ M nc . Here v stands for a multi-index 
of length n c . If F is sufficiently smooth, we may define 

B £ C 2 {X, X *), <fa)) := D 2 i?(0) (</>!, <fo) (4.2a) 

C £ £ 3 (X, X *), C(0i, 2 , 3 ) := ^ 3 i?(O)(0i, 2 , fa) (4.2b) 

for all (f>i G X. The nonlinearity -R : X — > X * may then be expanded as 

R(4>) = l -B{fa </>) + lc(0, 0, 0) + O(||0|| 4 ) (4.3) 

Let % : V C M nc — > X be a mapping that is as smooth as F and defined on a neigh- 
bourhood V of the origin in the coordinate space M nc with image H(V) = W^ c . Then % 
admits an expansion 

E ^h v z v + 0(\zf) (4.4) 

1<|i/|<3 

where z/ is a multi-index of length n c and /i„ £ I is an unknown coefficient. By the 
invariance of Wf oc we have 

%{z(t)) = u{t) VtGlR 

Differentiating both sides with respect to time leads to the homological equation 

A & *jU{z) + R(H(z)) = j{DU{z)z) (4.5) 

Substituting the expansions (4.1), (4.3) and (4.4) into (4.5) and equating coefficients of 
the corresponding powers of z, one recursively obtains the unknown coefficients h v and g v 
by solving linear operator equations of the form 

(A - A & *)(j) & * = </> * (4.6) 

where A G C and £ X * is given. If A o~(A) then (4.6) has a unique solution 
* £ D(A®*) for any given right-hand side. On the other hand, when A £ o~(A) a 
solution * of (4.6) need not exist for all right-hand sides ip & * . The following key lemma 
provides a condition for solvability that is useful in this situation. 

Lemma 33 (Fredholm solvability). Let A £ C \ {—a}. Suppose that L\ £ C(Y) defined 
in (3.7) is compact. Then A — A® : D(A®) C X — > X has closed range. In particular, 
(4.6) is solvable for * £ D(A G *) given ip & * £ X * if and only if = for all 

(p Q £ N(X-A*). 

Proof. From Corollary 15 in §3.1 we infer that TZ(X — A*) is closed. We first prove that this 
implies that 1Z(\ — A e ) is closed as well. Indeed, let (^ )neN be a sequence in 1Z(\ — A e ) 
such that V> — > V' ^ X e . Then there is a sequence (0 )neN m D(A e ) such that 

^ = (A - A & )^ = (A - A*)^ Vn £ N 
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where (2.14) was used in the second equality. Hence ?/> G 1Z(X — A*) for all n G N, so 
there exists G D(A*) such that (A - A*)(f>® = ip Q . Now 

A*<f> = -(A - A*)0 e + A0° = --0® + A0° G X 

so 0® G and (A - A & )4>® = ip & by (2.14). Hence V G TZ(X - A°). 

The second statement in the lemma is obtained from Banach's Closed Range Theorem 
[55, §VII.5] by which it follows that (4.6) has a solution if and only if annihilates 
N(X-A & ), i.e. if and only if 

(0 ,<0 *)=O V^e A^(A — A & ) 

To conclude the proof we show that iV(A - A G ) = N(X - A*). Indeed, N(X - A ) C 
N(X-A*) by virtue of (2.14). Conversely, suppose that 0° G N(X-A*). Then 0° G D(A*) 
and A*<p Q = A0° G X . Hence N(X - A Q ) D iV(A - A*) again by (2.14). □ 

4.2. The Andronov-Hopf critical normal form. In this case cr(A) contains a simple 
purely imaginary pair X\ t 2 = ±iu)Q with ujq > and no other eigenvalues on the imaginary 
axis. Let and 0® be complex eigenvectors of A and A* corresponding to Ai = iujQ and 
satisfying (0,0®) = 1. The restriction of (DDE) to the critical center manifold W,° is 
smoothly equivalent to the Poincare normal form 

z = iu z + g 21 z\z\ 2 + 0(\z\ 4 ) (4.7) 

where z is complex and the critical normal form coefficient g 2 i is unknown. Any point 
£ in the real two-dimensional center subspace Xq corresponding to Ai 2 may be uniquely 
expressed with respect to the set = {0, 0} by means of the smooth complex coordinate 
mapping 

The homological equation (4.5) presently becomes 

A & *jrl(z,z) + R{n(z,z))=j (D z n(z,z)z + DzH(z,z)f) 
with center manifold expansion 

-H{z,z) = z<j> + z4>+ ^ h jkz j z k + 0(\z\ 4 ) 

Note that since the image of H lies in the real space X, it follows that its coefficients 
satisfy hkj = hjk- The derivates z and z are given by (4.7) and its complex conjugate. 

Comparing coefficients of the quadratic terms z 2 and zz leads to two non-singular linear 
equations for jh 2 o and jh\\ with solutions 

jh 20 = -(A°*r i B(<t>,® 

jh n = (2ilo - A & *)~ l B(4>, 0) 

There are two equations corresponding to the cubic terms z 3 and z 2 z, the first of which 
is non-singular. The second one reads 

(iu I - A & *)jh 21 = C(0, 0,0) + £(0, h 20 ) + 25(0, h u ) - 2g 21 j<P (4.9) 

An application of Lemma 33 to (4.9) yields 

521 = ~(0 ,C(0,0,0) + £(0,M + 25(0,M> (4.10) 
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with /120 and h\\ implicitly given by (4.8). The cubic coefficient 521 determines the first 
Lyapunov coefficient l\ by the formula 

h = —Re 521 

It is well known [38] that in generic unfoldings of (4.7) l\ < implies a supercritical 
bifurcation of a limit cycle on the corresponding parameter-dependent locally invariant 
manifold, while l\ > implies a subcritical bifurcation of a limit cycle there. 

Remark 34. Notice that the vector (f>® satisfies A*<f>® = iujQcjJ 3 instead of A*(j) & = —icoo4> & 
which is used in the finite dimensional case. The reason for this is that the pairing (-, •) 
between X® and X & * is complex-linear in both arguments. Also, observe that the values 
of the multilinear form in (4.10) and (4.8) are elements of the dual space X & * of X®, i.e. 
they are (bounded) linear functionals. This has been taken into account in the numerical 
computations of §5. A similar remark is valid for the expressions in §4.3. 

4.3. The double Hopf critical normal form. In this case o~(A) contains two simple 
purely imaginary pairs 

A 1)4 = ±iui, A 2 ,3 = ±i^2 
with (jj\ t 2 > 0) an d no other eigenvalues on the imaginary axis. Let 4>\^ and 0® 2 be 
eigenvectors of A and A* , 

A(f>i = iu)i<fii, A(f>2 = iu2<t>2, A*4>f = iuji4>f, A*4>% = iw 2 ^f 

As in the finite-dimensional case, it is always possible to scale these vectors such that the 
'bi-orthogonality' relation 

{<t> h 4>f) =<% (l<i,i<2) 

is satisfied. In addition, we assume the non-resonance conditions 

kuji / Iuj 2 for all k, I e N with k + I < 5 (4.11) 

Then the restriction of (DDE) to the critical center manifold Wf oc is smoothly equivalent 
to the Poincare normal form 



( z\ = iuizi + 5210021 1 ^i 1 2 + gwiizi\z 2 \ 2 + #3200^1 1 2i 1 4 + g2inzi\zi\ 2 \z 2 \ 2 

+ 51022 21 1 2 2 1 4 + 0(\\(Z1,ZX, 2 2 ,22")|| 6 ) 
Z 2 = iuj 2 Z 2 + 5lll022|2l| 2 + 5002l22|22| 2 + 9221022|2l| 4 + 5ll2l2 2 1 2i | 2 1 Z 2 \ 2 
+ 50032 2 2 1 2 2 1 4 + O ( 1 1 (Zl , Zj, Z 2 , zt) \ | 6 ) 

where the constants gjkim are all complex [38, Ch. VIII]. Define 



(4.12) 



Pll 


P12 




52100 


5ion 


P21 


P22_ 




51110 


50021 



and assume that 

RepnRepi 2 Rep2iRep22 / 
As in shown in [38, Ch. VIII] the restriction of (4.12) to W!l is locally smoothly orbitally 
equivalent to 



z\ = iuj\Z\ + p n zi\zi\ 2 +p\2Z\\z2^ + iriZi\zi\ A + sizi|z 2 | 4 



+ 0(\\(Z 1 ,Z 1 ,Z2,Z2)\\ 6 ) 
P", 

+ 0(\\(Z 1 ,Z 1 ,Z2,Z2)\\ 6 ) 



Z 2 = iu) 2 Z 2 +P2\Z2\z\\ 2 +P2222|22| 2 + S 2 Z2\zi\ 4 + ^2-22 1 22 1 4 



(4.13) 
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Here pij and Sj are complex while t% are real, for 1 < i,j < 2. The real parts of Si are 
given by 



Re 51121 _ Re 50032 _ Re 53200 Re 50021 
Re 5mo Re 50021 Re52iooRe5mo 

Refill n Re 53200 Re 52iooRe 50032 



Resi = Re 51022 + Re 51011 X 

and 

Re s 2 = Re 52210 + Re 51110 x 

Ke 51011 Re 52100 Re gioiiRe 50021. 
The real constants Ti are of secondary importance in the bifurcation analysis of a generic 
two-parameter unfolding of (4.13) and so we omit expressions for these. They can be 
extracted from the proof of [38, Lemma 8.14]. 

The double Hopf bifurcation is a complicated bifurcation, both from a computational 
as well as a conceptual viewpoint. An unfolding of (4.13) is best analysed by rewriting it 
in polar coordinates. The sixth-order terms may not be truncated, since they may affect 
the qualitative dynamics. Depending on the sign of 

RepiiRep 2 2 = Re 52iooRe 50021 

this bifurcation exhibits either 'simple' or 'difficult' dynamics, see [38, §8.6.2]. Assuming 
generic dependence on parameters, one may encounter invariant tori, chaotic dynamics, 
Neimark-S acker bifurcations of cycles and Shilnikov homoclinic orbits. Note that, al- 
though computations up to and including fifth order are required to determine all critical 
coefficients, computations up to and including third order suffice to distinguish between 
'simple' and 'difficult' cases. 

The critical normal form coefficients may be obtained using a procedure similar to the 
Hopf case discussed in §4.2. We omit the details and only present the results, noting that 
the center manifold now has the formal expansion 

H(Z1, Z\, Z 2 ,Z 2 ) = £101 + ^101 + 2 2 2 + #202 + TmTi - ] hjkl m z{z'{Z 2 Z 2 n 

j+k+l+m>2 J 

At the second order in the corresponding homological equation, we find 

jhnoo = -(^T^i,^) 

jh 20 oo = (2iu 1 -A e *y 1 B(4> 1 ,(l )1 ) 

jh 1010 = [i{u)i + u 2 )-A Q *]- 1 B(<i> 1 ,<l) 2 ) 

j/iiooi = [*(wi - u) 2 ) - A®*] -1 B((j>\, 2 ) 

jh 0020 = (2iuj 2 - A & *Y 1 B((I) 2 , 2 ) 

All operators in the right-hand side of the above equations are invertible due to the 
assumptions (4.11) on the critical eigenvalues. 

Further, one obtains the following equations for hjkim with j + k + I + m = 3: 

j/isooo = (3zwi - A®*)- l [C{(t> u 0i, 0i) + 35(^000, 00] 

jh 2010 = [i(2ui + u 2 ) - A *] ~ l [C (0i , 01 , 02) + B(h 200 o, <h) + 2B(frioio, 0i)] 

jh 20 oi = [i(2wi -u 2 ) - A *] _1 [C(0i, 0i, 2 ) + B(h 200 o, 2 ) + 2B(/iiooi,0i)] 

jh W20 = [i(wi + 2w 2 ) - A G *]- 1 [C(0i, 2; 2 ) + B(h 0020 , 0i) + 2B(h mo , 2 )] 

jhum = [i(ui - 2uj 2 ) - A 0Jr ]- 1 [C(0i,0 2 ,0 2 ) + £(W),0i) + 2B{h 100 iA 2 )] 

jhooso = (3io; 2 - A *)" 1 [C(0 2 , 02, 02) + 35(^0020, 02)] 
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The cubic coefficients in the normal form (4.12) come from the Fredholm solvability con- 
ditions and are given by 

52100 = 'C^i'^i'^i) + ^( /l 20oo,<M + 2 J B(/inoo,<M) 

51011 = (0f,C(0l,0 2 ,0 2 ) + ^(/lioio, 2 ) +-8(^1001,02) + ^(/lOOll, 0l)) 
51110 = (4>2 ,C(01,01,02) +-B(/lnoo, </>2) + -6(^1010, 0l) + #(^1001, 01 )) 

50021 = 2 1 (02 9 ,C(0 2 ,02,02) + -B(/lOO2O,0 2 ) + 2 J B(/l O Oll,02)) 

Similarly, one can compute all remaining coefficients in (4.12) by proceeding to orders four 
and five. The resulting (lengthy) formulas are omitted. For the finite-dimensional case 
these can be found in [37]. 

4.4. Evaluation of normal form coefficients. The computability of the normal form 
coefficients derived in the previous subsections depends on the possibility to evaluate 
the dual pairing (0®,0®*), where 0® G X® is some eigenvector of A* corresponding 
to a simple eigenvalue A G <j{A) and * G X Q *. Moreover, the coefficients h u , with 
v a certain multi-index, can only be computed once a representation for the resolvent 
R(X,A Q *) is known, where A 6 p(A). At first sight this seems to be a difficult task, since 
X® = Y* x Z+QO, h];Y*) and hence 

X e * = Y** x [L^h^Y*)]* 

see §2.2, and, as remarked there, [-^([O, h]; Y*)]* / L°°([—h,0];Y**). Moreover, a repre- 
sentation of the second dual space Y** is generally unknown, e.g. when Y = C(f2) as for 
(NF). 

Remark 35. In §5 it will turn out that the second derivative B in (4.2a) vanishes due to a 
symmetry in (NF) for the particular modelling functions chosen. In the present subsection 
we deliberately do not exploit this information in order to illustrate a general principle. 

In this subsection we offer a way around these complications that works for equations of 
the type (DDE). We first deal with the problem of determining R(X, A 13 *). From Lemma 
4 in §2.2 it follows that the second and third derivatives defined in (4.2a) and (4.2b), as 
well as all derivatives of higher order, map into the closed subspace Y x {0} of X®* . By 
inspection of the expressions for the coefficients h u in §§4.2 and 4.3 one sees that it is 
sufficient to obtain a representation of the action of R(X, A®*) on this space. 

Lemma 36. Suppose that A G p(A). For each y G Y the function ip = E\ A(X)~ 1 y is 
the unique solution in C 1 ([— h, 0]; Y) of the system 

f AV(0) - DF(Q)iJ; = y 

Moreover, ■0®* = jtp is the unique solution in D{A®*) of (A — A Q *)ip & * = (y, 0). 

Proof. We return to the setting of Proposition 14 in §3.1 with L = DG(0). Since A G p(A) 
it follows that A(A) _1 exists. We start by showing that ip = E\ ® A(A) _1 y solves (4.14). 
Explicitly, 

^{d) = e xe /\{X)- l y V6e[-h,Q] 
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so clearly ip G C 1 ([— h, 0]; Y) and ip satisfies the second equation in (4.14). Prom (2.5) we 
recall that PP(0)V> = -atp(0) + DG(0)ip. Therefore, 

Xip{0) - DF(0)^j = (A + a)V»(0) - £>G(0)V> 

= (A + a)A(A)- x y - P>G(0)(e A A(A)- X y) 

= (A + a)A(A)- 1 y-L A (A(A)- 1 y) 

= A(A)A(A)- 1 y = y 

2.2 implies that jij) G D(A & * 



Lemma 9 in 
and 



(\-A°*)j^ = \ 



But a(A & *) = a (A) so -0®* = jif> is the unique solution of (A — A *)-0 * 
Consequently, ip itself is the unique solution in C l ([—h, 0];Y) of (4.14). 



where j is the embedding defined in (2.8), 

(y,o) 



(y,o). 

□ 



" ^(0) " 




"£>F(0)?/>~ 









The above lemma takes care of one of the two problems sketched above. Now suppose 
that A G cr(A) \ {—a} is a simple eigenvalue with eigenvector <p G -D(j4). (Note that A is 
isolated in cr(A) by Corollary 18 in §3.1.) Let 4> e G D(A*) be a corresponding eigenvector 
of A*. Without loss of generality we may assume that {(f>,4>®) = 1. Let P and P * be 
the associated spectral projections on X and X *, respectively. We set out to evaluate 
(0°,(/) *) where <fi & * = (y, 0) G Y x {0} C X * is given, but <^ is unknown. Since the 
range of P * is spanned by j<f> we have P©*<^©* = Kjcj) for a certain « G C. In fact, by 
(2.8) it follows that 

((j> Q ,(j) & *) = (P </> ,</> *) = (<j) & ,P & *<p Q *) = k(^ g J(/)} = k (4.15) 

so k is to be determined. This may be done as follows. From the Cauchy integral repre- 
sentation [45, §5.8] for P * we infer that 

1 



0*^0* 



2m 



dc x 



R{z,A°*)</)°*dz = Kjc 



(4.16) 



where C\ is any open disk centered at A such that C\$ C p(A) where Ca,o : = C\ \ {A} 
and dC\ is its boundary. Since (fi G * G y x {0} the integrand in (4.16) may be calculated 
using Lemma 36. Specifically, for z G dC\ we have 

R(z, A Q *)(f) Q * = j(e z A(z)- X y) = Ai ' r> ' 



y 



(4.17) 



_ £^«>A(z)- 1 y 

Since j<f> = (j>(0) we may restrict our attention to the first component to infer that 

^-l A(z)- 1 ydz = «0(O) 
27 ™ Jdc x 

We note that the integral is l"-valued and (4.17) is an identity in Y. The integrand may 
be evaluated using the results of §3.4. Indeed, for each z G dC\ it is necessary to solve a 
system of the type (3.34), but with S z cp replaced by y. 

For this purpose we may apply Theorem 30 as follows. Let us assume that A is a root 
of the characteristic equation (3.32) on the imaginary axis, A S, the roots ±/?j(A) of the 
polynomial (3.26) are all distinct and (3.29) holds. Suppose it has also been verified that 
the matrix T(A) is invertible. By choosing the radius of C\ sufficiently small, we guarantee 
that for every z G C\ it holds that z ^ —a, z G" 5, the roots ±pi(z) are all distinct, (3.29) 
is satisfied (with z instead of A) and T(z) is invertible. Furthermore, in this way we may 
also ensure that z G p{A) for every z G C\,o- In particular, all the maps 

C x 3 z t-> ± Pi (z) G C (i = l,...,N) 
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are analytic. By (3.39) and (3.40) (with y in place of h z ) this implies that 

C x Bz^ f z (x) G C 2N 
is analytic for all x G fi. Hence by (3.42) (with y instead of /i 2 ) we have, for every x G f2, 



9C7 A 



2/0*0 
0C\ 2 + « 



A' 



i=l 



i=1 L Jac\ Jdc\ 



(4.18) 



where To i2 is as in (3.41). Since ro j2 involves S(z) 1 , the maps 

CA3^[r ,d±ieC (i = i,...,JV) 

cannot be expected to be analytic and (4.18) may not be reduced further. 

In summary, (4.18) provides a way to evaluate (4.17) by numerical integration. It 
suffices to parametrise dC\ and apply a quadrature rule to compute the C- valued contour 
integrals 



\To,z]±idz (i = l,...,N) 



which are independent of x G Q. We then verify that i times (4.18) and 0(0) are indeed 
proportional to each other as functions of x G O. The value of (0®,0 *) in (4.15) then 
equals the corresponding proportionality constant n. 



5. Numerical calculations 

In §3.3 we derived a characteristic equation for problem (3.4) under the assumption 
that J is a finite linear combination of exponentials. The main result was formulated 
in Theorem 28. Subsequently, in §3.4 we obtained a closed expression for the associated 
resolvent operator. In the present section we apply these findings together with the theory 
from §4 to a concrete example. For reasons that will become apparent later, we assume 
that the connectivity function has a bi-exponential form, 

J(x,r) = cie-^ |x - r| + c 2 e-^ x - T \ Vx.refl (5.1) 

and we choose the activation function S as in [27], 



S{V) = - ^t-^ VT/G 

1 + e~ rV 2 



Since S(0) = it follows that (NF) admits the trivial steady state V = on which we will 
focus from now on. Here we have S'(0) = \ and hence Cj = j£j for i = 1, 2. 

Let us continue by expressing the characteristic equation for this example and discussing 
a naive approach for finding its roots. Thereafter we compare these results with a more 
traditional approach which discretises the spatial domain £1. Such a discretisation can 
be studied using techniques and software that are already available. We conclude with a 
normal form analysis of a Hopf bifurcation and a double Hopf bifurcation to illustrate the 
potential of the results from §4. 
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5.1. Spectral calculations. In order to apply Theorem 28, we start by considering the 
characteristic polynomial V from (3.26), which presently takes the form 

V( P ) = eXT ° {X 2 + a \ p 2 (A + M!) 2 )(P 2 - (A + u 2 f) 

+ ci(A + ^i)(p 2 - (A + /i 2 ) 2 ) + c 2 (A + p 2 )(p 2 - (A + mi) 2 ) 

This is a second order polynomial in p 2 . We apply a Newton algorithm to the mapping 
A i y det5(A) to find the solutions of the characteristic equation (3.32). At each root A 
we need to verify that A S, the numbers ipi^A) are all distinct and (3.29) is satisfied. 
(Note that both of these are open conditions.) Passing this test we may conclude that A 
is indeed an eigenvalue. 

5.2. Discretisation. 

Derivation [27] . An approximate solution to the neural field equation can be obtained by 
discretising (NF). This reduces the state space from C([-h,0];Y) to C([-h, 0]; R m+1 ) 
for some m G N. Hence the theory of 'classical' DDE can be applied to analyse the 
approximate system. 

We heuristically derive the discretised system following [27] but we make a few minor 
corrections. Consider the original equation (1.1): 

dV f Xi 

— (t, x) = -aV(t, x) + 22 J ( x ' r ) S ( v ( t ~ T ( x > r )> r )) dr 
dt i=i Jx *-i 

for some partition — 1 = xq < x\ < . . . < x m = 1. We approximate every single integral 
with a two-point trapezoid rule evaluated at the end points of the integration interval, 

dV m x — x 

— (t,x) « -aV(t,x) + ^2 - l ~ [J(x,Xi-i)S(V(t - T(x,Xi-i),Xi-i)) 

i=i 

+J(x, Xi)S(V(t - t(x, Xi),Xi))] 
By writing Vj(t) = V(t, Xj), we obtain for j = 0, 1, . . . , m, 

-^(t) = -aVj(t) + Y, Xt [J(xj, Xi_i)5(^-i(t - rfo,^))) 

i=l 

+ J(Xj,Xi)S(Vi(t - T(Xj,Xi)))] 

As in (3.14) we take r(x,r) = tq + \x — r\. Also, with some abuse of notation we write 
J(\x — r|) for J(x,r), since the dependence in the right-hand side of (5.1) on x,r is only 
via \x — r\. By restriction to an equidistant mesh of size 5 = x% — = ^, we obtain 

dV- 9 m 1 

~^{t) = -dVj(t) + - E 2 W " j ~ 1 D 5 (^-i(* -r -6\i-j-l\)) 

+j(S\i-j\)S(Vi(t-T -6\i-j\))] 

Defining 

\\ ifiG{0,m} 

Wi = { r 

[1 if i e {l,2,...,m- 1} 
enables us to telescope the summation, arriving at 
dV- 2 m 

-£(t) = -aV.it) + WW ~ 3\) S ( V ^ ~ to ~ S \* ~ iD) (DNF) 
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Figure 1 . Comparison between spectra of the discretised system for m = 
20, 50, 100 and roots of the characteristic equation. The four black arrows 
indicate four roots which are not in the spectrum and the grey arrows point 
out distinct values which are not found. See text for a more elaborate 
description of these points, a = 1, ro = 1, c\ = —5, C2 = 2, [i\ = 2, = 0. 

for all j = 0, 1, . . . , m. We refer to (DNF) as the discretisation of (NF) or (1.1). Note that 
(DNF) indeed is a classical DDE, albeit with many delays, which may be implemented in 
MATLAB to perform forward-time simulations using the dde23 scheme. In particular, the 
software package DDE-BIFT00L [2 ' } allows us to determine the spectrum of the discretised 
system. At the end of this section we consider two examples in which we use both our 
analytic results and these numerical tools to study critical points in neural fields. 

Convergence of discretisation. In order to 'validate' the above discretisation procedure, 
we generate discretisations with different resolutions and compare their spectra with the 
spectral values obtained by using the methods from §3.3. This is illustrated in Figure 1. 
The black arrows indicate four roots of the characteristic equation which do not satisfy all 
conditions stated in Theorem 28: these points do satisfy (3.32), but also pi(A) = /02(A). 
Therefore they are to be rejected as eigenvalues. Furthermore, the grey arrows indicate 
eigenvalues which are not found using the algorithm of §5.1. These eigenvalues all lie in 
the 'accumulation region' near the point —a. 

We proceed by studying this accumulation of eigenvalues more closely. For A f — 1 along 
the real axis the absolute value of detS'(A) is plotted in Figure 2 on log-log scale. Every 
downward peak corresponds to a root of the characteristic equation. For A close to —a the 
numerical accuracy drops causing the peaks to be less pronounced. This shows that, while 
these are not located by the root finder algorithm, the characteristic equation does have 
accumulating roots near the essential spectrum {—a} as is suggested by the spectrum of 
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1+1 



Figure 2. Detail of the accumulation of eigenvalues along the real line near 
A = — 1: | det(S*) | is plotted near the accumulation point and downward 
peaks correspond to roots of the characteristic equation. 

Table 1 . Parameters corresponding to Hopf bifurcation 



parameter 


a 


TO 


r ci 


C2 


Mi 


M2 


value 


1.0 


1.0 


4.220214885988226 3.0 


-5.5 


0.5 


1.0 



the discretisation, cf. Figure 1. However, due to both the high frequency oscillations and 
numerical errors, the Newton algorithm is unable to locate these roots. 

Finally we observe that spectra corresponding to finer meshes converge to the analytic 
spectrum. However, it appears that for increasing resolutions DDE-BIFT00L focuses on 
roots near —a instead of eigenvalues located further away. This is clearly seen when 
m = 100, in which case no spectral values A are found with Re A > —1.2. For that reason 
we have chosen m = 50 in the following examples. 

5.3. Hopf bifurcation. Rhythms and oscillations are important features of nervous tissue 
that could be studied with neural field models. For that reason Hopf bifurcations play a 
key role in the analysis of neural field equations. In this subsection we study a concrete 
example of a Hopf bifurcation, both analytically and numerically. We also compare the 
results of both methods. 

Initially we focus on a connectivity of the 'inverted wizard hat '-type. Similarly as in 
we consider the steepness parameter r of the activation function as bifurcation parameter. 

Spectrum. The characteristic equation (3.32) is used to determine the point spectrum for 
a range of parameters. For the values shown in Table 1 there exists a purely imaginary 
pair of simple eigenvalues. The corresponding spectrum is displayed in Figure 3. 
The figure also shows the spectrum as calculated by DDE-BIFT00L for a discretisation 
of m = 50 intervals. From the graph it is obvious that the solution algorithm of §5.1 
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O Discretisation 
x Analytical 



Figure 3. Spectrum at a Hopf bifurcation. Comparison between analytic 
approach and discretised system (m = 50). 



is unable to locate eigenvalues near the accumulation point —a = — 1. We discuss this 
phenomenon below. Apart from that, the numerical approximation seems to be very close 
to the analytic solution. Only in the far left half-plane an error can be observed. From a 
dynamical point of view such an error is of course rather innocuous. 

First Lyapunov coefficient. In order to determine analytically the type of Hopf bifurcation 
(i.e. sub- or supercritical), the first Lyapunov coefficient has to be determined. Before 
the result of §4.2 can be applied, the eigenfunction corresponding to the eigenvalues at 
criticality has to be determined. Application of Theorem 28 yields 



4>(t,x) 



[ 7l (ePi* + e - 



P2X\ 



Vt e [-h,0], Vxen 



with 



pi = 0.321607348361597 - 0.880461478656249i 
p 2 = 0.110838003673357 - 2.312123026384049z 

71 = -0.191821747840362 - 0.172140605861736i 

72 = -0.080160108888561 

corresponding to A = 1.6440031020468931 (Note that in the present example with to as 
in Table 2 and Q = [—1, 1] the delay interval equals [—h, 0] = [—3, 0].) Since the activation 
function S is odd, its second derivative vanishes and the critical normal form coefficient 
</2i in (4.10) significantly simplifies to 

521 = ^{(f> & ,D 3 R(0)(4>,(f>,^)) 

The pairing is expressed as a contour integral around A which we evaluate numerically, 
see §4.4. We find (721 ~ —0.326 + 0.0389? and hence the first Lyapunov coefficient is 
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Figure 4. Forward time simulation of discretised system (m = 50) for 
r = 6 beyond a Hopf bifurcation. A long transient is observed before the 
solution approaches the limit cycle. 



Table 2. Parameters corresponding to double Hopf bifurcation. 



parameter 


a 


TO 


r ci 


C2 


Mi 


M2 


value 


1.0 


1.0 


4.828749714457348 3.0 


-5.5 


0.0 


0.999592391420082 



Table 3. Values of A and corresponding p 



A 




2.030930500644927i 


0. 454550410967142-1. 057267648955222i 
0.054136932895367-3.495632804443535i 


1.299147304907829i 


1.075429529957343-0.717519976488838? 
1. 128716151852882-2. 306528729845143i 



li k, —0.198. The negative sign of l\ indicates a supercritical Hopf bifurcation. Therefore 
stable periodic solutions are expected to emerge from the bifurcating steady state. 

Simulations. We choose r = 6 which is beyond the critical parameter value of Table 1. 
For the initial condition V(t,x) = e = 0.01 for t £ [—h,0] with h = 3 and x £ [—1,1] 
the simulation result is shown in Figure 4. After a transient time, the system approaches 
its stable periodic attractor. The convergence to stable periodic motion is consistent with 
the sign of the first Lyapunov coefficient. Furthermore, both the shape and period of this 
attractor match with the eigenfunction and eigenvalue respectively. 

5.4. Double Hopf bifurcation. The spectrum at the Hopf point studied in §5.3 con- 
sists mainly of complex pairs of eigenvalues. Therefore it is to be expected that system 
parameters can be tuned such that a second pair of complex eigenvalues arrives at the 
imaginary axis, giving rise to a double Hopf bifurcation. In this subsection we show that 
this is indeed possible and we study this bifurcation both analytically and numerically. 

Spectrum. Parameters for which the system has two complex pairs of eigenvalues on the 
imaginary axis and no eigenvalues in the positive right half-plane are identified, see Table 
2. The corresponding spectrum is depicted in Figure 5 while Table 3 lists the values of 
A at this critical point. As with the regular Hopf bifurcation, we observe that the root 
finding algorithm misses most eigenvalues near the essential spectrum at —a = —1. 
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Figure 5. Spectrum at a double Hopf bifurcation. Comparison between 
analytic approach and discretised system (m = 50). 
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Next we compute the eigenfunctions corresponding to the critical eigenvalues. For this 
purpose Theorem 28 may be applied using the data in Table 3. Modulus and argument of 
both eigenfunctions are depicted in Figure 6. 

Normal form coefficients. With this information available, the normal form coefficients of 
the double Hopf bifurcation may be evaluated. The coefficients 92100, flloiij 91110, and 
90021 from §4.3 are found as in §5.3, which results in the matrix 

-8.822 -3.367" 
-13.79 -1.310 

Since P11P22 > 0, we conclude that this double Hopf bifurcation is of the 'simple' type, 
see [38, §8.6.21. Defining 9 := ^ » 2.57 and 5 := ^ » 1.56, we find that 96 > 1 and 
therefore this 'simple' bifurcation has sub-type I [38]. The key feature of this sub-type is 
the presence of a regime in parameter space for which two distinct stable periodic solutions 
exist. 

Simulations. The parameters are adjusted slightly, such that both pairs of complex eigen- 
values have a positive real part. More specifically we choose r = 6 and fi2 = 1 while 
keeping other parameters as in Table 2. For the following two initial conditions 

V(t,x)=ex (5.2a) 
V(t,x)=e (5.2b) 

with s = 0.01, the discretised system (m = 50) is integrated forwards in time, see Figure 
7. 

For the chosen parameters the system has two stable periodic attractors and the asymp- 
totic behaviour is determined by their initial conditions. This result is consistent with the 
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Figure 7. Bi-stability near double Hopf bifurcation: for r = 6 and [1% = 1 
the time evolution is shown for different initial conditions (m = 50). Top 
and bottom diagrams correspond to (5.2a) and (5.2b) respectively. 
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predictions of the normal form computation. Furthermore we observe, since the system is 
close to the double Hopf bifurcation, that both the shape and period of the periodic so- 
lutions are fairly well approximated by the critical eigenfunctions (c.f. Figure 6). Indeed, 
the moduli of the eigenfunctions correspond to the amplitude of the asymptotic solution. 
For either of the solutions the extrema are located near the borders of the domain. The 
antiphasic solutions in the upper panel are indicated by the jump of size n in the argument 
of the first eigenfunction, see Figure 6 on the left. 



6. Discussion 

We have demonstrated that neural field equations with transmission delay fit well into 
the sun-star framework for delay equations. As a consequence, standard results from 
dynamical systems theory, such as the principle of linearised (in)stability, center manifold 
reduction and normal form computation, are readily available. These, in turn, open the 
possibility for a systematic study of codimension one and two local bifurcations, w.r.t. 
parameters in the connectivity and activation functions. This facilitates an understanding 
of the effect of parameters in terms of biological quantities. 

In §5 we analysed the dynamics of a one population model with the inverted wizard 
hat as connectivity function. The choice of an inverted Mexican hat is biologically more 
plausible, as pyramidal cells are surrounded by a cloud of interneurons, while long range 
connections are by and large excitatory. We have chosen the inverted wizard hat mainly 
for mathematical convenience. Indeed, in §3.3 an analytic formula for the location of the 
eigenvalues was derived. It is well known that the combination of an inverted Mexican 
hat connectivity with a transmission delay leads to dynamic instabilities [6, 34]. In [12] 
Turing instabilities were shown to occur for the inverted wizard hat connectivity. 

We have seen that the stationary spatially homogeneous state destabilises upon in- 
creasing the steepness (gain) of the activation function. This is in line with other findings 
indicating that the activation function strongly influences dynamical behaviour, see for 
instance [26, 13]. 

It is mathematically challenging to consider neural field equations on unbounded spatial 
domains, leading to infinite delays, although such is of less importance from a biological 
viewpoint. Our main goal for the near future is to develop tools for numerical bifurcation 
analysis for the class of equations studied in this paper. Normal form computation is a 
first prerequisite for this task. Hence we are on our way. 



Appendix A. Proof of Proposition 26 

We use the same notation as in Lemma 23 and its proof. We recall that the vector 
Z = [Co> Cl> • • • j CiV-ij 1] is chosen such that the vector f3 = [flo, (3\, . . . , f3jsi], whose elements 
are the coefficients of the characteristic polynomial V, is given by f3 = M T Z. Introducing 
r : = [1, p 2 , p 4 , . . . , p 2N ] we see that 

V{p) = r T M T Z (A.l) 

First we determine the vector Z, thereafter we split M, and we conclude the proof by 
determining how Z acts on each factor in this decomposition. 
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Although Z can be obtained by applying the inverse of the Vandermonde matrix W, 
we will proceed in a different manner. We start by rewriting (3.22) as 





"1 k 2 


1 •• 




1 k\ 


k\ 




1 k 2 


k 4 
N ■ ■ 







. . 


define P m : 


= [Pl,P2 



k\ N ~ 2 0" 

kl N ~ 2 



k 2 N N - 2 o 





" Co " 




~-k 2N ~ 




Ci 




k 2N 




Cjv-i 




-k 2N 




1 




1 



(A.2) 



\Pm\ = Y^JiLiPi equal to the number of l's in P m . Using Gaussian elimination the solution 
of (A.2) is found to be 



V i \N-o v h 2pi h 2p2 h 2pN 

1 ' 2^\P N \=N-0 K l K 2 ' ' ' N 



I 1\W-lv^ ,2pi,2p 2 r,2 Pjv 

I L ) Lu\P N \=N-\ ft l K 2 ■ ■ ■ K N 



( I^V h 2pi h 2p2 h 2pN 

\ L ) Z_JP W |=1 ft l K 2 ■■■ K t 



1 



From the proof of Lemma 23 we recall the decomposition 

M T = e Xr °(X + a)I + 2E 



(A.3) 



where / is the (N + 1) x (N + 1) identity matrix and H was defined in the proof of Lemma 
23. Expanding the bilinear forms in the matrix S and moving the summation in front of 
the matrix yields 



N 



n — ^ Ciki^i 



1 kj 

1 







„2(JV-1)- 
2(N-2) 



(A.4) 



Now substitute (A.3) with (A.4) into (A.l) to obtain 

N 



T(p) = r T e Ar °(A + a) 1 + 2^2^3 

i=l 

= e Ar °(A + a) [l p 2 p 4 ... p 2N ]Z + 2r T Y,Cik l ~ i Z 
We observe that on the one hand, 

e AT0 (A + a)[l p 2 p 4 ... p 2N ] Z = e XT0 (\ + a)Y[{p 2 -k. 



N 



i=l 



N 



(A.5) 



i=i 
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while on the other hand, 

N 



N 



5>^[i p 4 ... p™] 



i=l 



~ L ) 2^\P N _ 1 \=N-1 ^1 ft 2 • • • K i+1 ■ ■ ■ K N 

-,\N-2ST ,2pi,2p 2 i?Vi-\ -L&i l?PN-l 

~ L ) Z^|P JV _ 1 |=AT-2 K \ ft 2 ■■■ K i-\ • • • 



1 





Hence by (A. 5) it follows that 

N N N 

V(p) = e AT0 (A + q) l[(p 2 - k]) + 2 ^ c^i [](p 2 - fc; 

i=l 4=1 i=i 



which is equivalent to (3.26), in the sense that the two polynomials have identical roots. 
Hence the proof is complete. 
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